Harmonic oscillator force between heavy quarks 
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A renormalization group procedure for eflective particles is applied to quantum chromodynamics 
of one flavor of quarks with large mass m in order to calculate light-front Hamiltonians for heavy 
quarkonia, Hx, using perturbative expansion in the coupling constant a\. A is the renormaliza- 
tion group parameter with the interpretation of an inverse of the spatial size of the color charge 
distribution in the effective quarks and gluons. The eigenvalue equation for H\ couples quark-anti- 
quark states with sectors of a larger number of constituents. The coupling to states with more 
than one effective gluon, and interactions in the quark-anti-quark-gluon sector, are removed at the 
price of introducing an ansatz for the gluon mass, /i^. The simplified equation is used to evaluate 
a new Hamiltonian of order a\ that acts only in the effective quark-anti-quark sector and in the 
non-relativistic limit turns out to contain the Coulomb term with Breit-Fermi corrections and spin- 
independent harmonic oscillator term with frequency uj = [(4/3)(QA/7r)]^''^ A(A/m)'^(7r/f f 52)^'^^. 
The latter originates from the hole excavated in the overlapping quark self-interaction gluon clouds 
by the exchange of effective gluons between the quarks. The new term is largely independent of the 
details of /i^ and in principle can fit into the ball park of phenomenology. The first approximation 
can be improved by including more terms in H\ and solving the eigenvalue equations numerically. 

PACS numbers: 12.38.-t,12.39.-x,12.90.-|-b,11.15.-q 



I. INTRODUCTION 

The purpose of this article is to describe a proce- 
dure that starts from a quantum chromodynamics (QCD) 
with only one flavor of massive quarks and produces the 
Schrodinger equation for heavy quarkonia in a single for- 
mulation of the theory. Only the first approximation for 
the final Hamiltonian is evaluated. In this simplest ver- 
sion, the procedure involves a guess for the gluon mass 
term. But the guess appears to have little influence on 
the result. The procedure itself is not limited to the sim- 
plest version and the gluon mass ansatz can be tested in 
future in refined calculations and phenomenology. The 
procedure is relativistic and can be used for quarkonia in 
arbitrary motion, which is a pre-requisite for application 
in high-energy processes. Chiral symmetry is explicitly 
broken in the case of heavy quarks and the issue of the 
spontaneous breaking of the symmetry is ignored. 

The approach described here stems from the similar- 
ity renormalization group procedure for Hamiltonians , 
which has been applied to QCD [3| using the light-front 
(LF) form of dynamics |3j|- A new ingredient is the boost- 
invariant creation and annihilation operator calculus for 
effective quarks and gluons (see below). Otherwise, the 
LF approach is known for a long time, mainly as a can- 
didate for connecting two qualitatively different models 
of hadrons: the parton model in the infinite momentum 
frame (IMF) and the constituent model in the rest 

frame of a hadron fol . Many contributions in that area 
[3 have followed the seminal work on exclusive processes 
[3. Through the boost invariance and precisely defined 
notion of effective constituents, the approach described 
here aims at providing a bridge between the two models 
of hadrons in a single theoretical framework. The case of 
a heavy quarkonium is chosen here as the simplest one 



to begin with and test the method. 

In the LF dynamics, the evolution of states is traced 
along the direction of a light-like four- vector n^, for which 
n? = 0. With the conventional choice of n = (1, 0, 0, — 1), 
the variable nx = -\- = x~^ plays the role of time 
while — p'^ ~ is the Hamiltonian. In order to 
define it for bare particles in QCD, one has to choose a 
gauge. No serious alternative exists to nA = A'^ = 0. 
But the equation D^F^^^ = implies a constraint that 
is analogous to the Gauss law and forces the Fourier com- 
ponents of A~ to contain inverse powers of the kinemati- 
cal momentum fc+ . Since k'^ ranges from to infinity, the 
inverse powers of fc+ produce singularities in the region 
around zero. 

One can impose a lower bound on fc"*" , such as A:+ > , 
to regulate the theory 0, 0, 0, ^ . The parameter (5+ 
becomes a smallest unit of momentum that any particle, 
physical or virtual, can carry in such discretized theories. 
But the fixed unit breaks the boost invariance required 
for connection between the IMF and the rest frame of 
any hadron. Namely, when some physically relevant P+ 
is made large, the smallest allowed x = /P~^ becomes 
small. In the IMF P+ 00 and the same small-x di- 
vergences re-appear despite the presence of S~^. The key 
singularity is related to the dx/x distribution of gluons 
in the parton model and seems to require a dynamical 
mechanism to remove. One cannot just vary together 
with P+ because boosts cannot change the cutoff in a 
quantum theory constructed ab initio |l2l |. 

There exists a possibility that the small-A:^ singular- 
ities are related to the properties of the vacuum state. 
The sum rules for h eavy quarkonia |l3l | include quark and 
gluon condensates that may participate in the dy- 
namics in the small-fc^ region 15, 16, 17, 1&, 19, 2,^. 
In the QCD picture with such nontrivial ground state 
I21L I23I , and relativistic bound-state excitations of this 
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state in the form of Tr-mesons, one can hardly hope to re- 
solve the small-/c^ singularity easily. The situation sim- 
plifies a lot in the case of quarks with mass m 3> ^qcd ■ 
The small- a: singularity and a non-perturbative bind- 
ing mechanism for quarks and gluons can interplay with 
each other without interference from the vacuum. In the 
case of light quarks, a similar interplay may be at work 
and contribute to the saturation mechanism of partons 
[23, 0, EE 113 J but that issue is not addressed here. 

The gluon mass ansatz is introduced to represent ef- 
fects of the non-abelian interactions. The ansatz is in- 
serted at the level of solving the eigenvalue equation 
for the Hamiltonian Hx, where A is the renormalization 
group (RG) parameter. The procedure of introducing the 
mass ansatz is similar to the one proposed in ^ and later 
discussed in simplified matrix models ,22j »28] . New ele- 
ments are the limitation of power counting to the relative- 
motion variables, the exact boost invariance, removal of 
small-x divergences through the mass ansatz as a func- 
tion of the relative momenta, and no need for ad hoc 
potentials in the first approximation. All these features 
will be described in detail later. 

In brief, H\ = T\ + V\, where T\ is the kinetic en- 
ergy operator made of all terms that are bilinear in the 
creation and annihilation operators for the effective par- 
ticles, and V\ represents all other terms. The parameter 
A defines the width of momentum-space form factors in 
V\. For some value of A = Aq ~ 1 GeV, one can freely 
add to a term of the form [1 — (aAo/as)"]M^j where 
n > 2. The number 2 ensures that corrections to the 
first approximation occur first in the fourth-order of ex- 
pansion of Hxo in powers of go = gxo, = 3o/(47r). 
This is the lowest order at which a perturbative shift in 
the gluon energy in any state can influence the contri- 
bution of that state to the dynamics of any other state 
in perturbation theory, keeping intact the QED-like small 
coupling expansion scheme with a Coulomb potential. 
stands for the mass term that one assigns to the effective 
gluons of the transverse size l/Ao- The interpretation of 
1/A as the spatial size of the color charge distribution 
in the corresponding effective particles is based on the 
feature mentioned above that Vx contains vertex form 
factors of width A in momentum space. The mass ansatz 
contributes to the invariant masses through /i^/x, where 
X is the fraction of the longitudinal momentum carried 
by the gluon. as denotes the large, relativistic value of 
the coupling constant in QCD at the scale Aq with a true 
^QCD in this scheme. The term with /i^ vanishes for 
ao = ag. Nevertheless, only /i^ counts when the ratio 
ao/cts is small. Thus, in the weak coupling expansion, 

Hx, = Tx, + + [Vxo - (ao/a,) V2] (1) 

and the term [Vxg — (ao/as)"^^] is treated as a source of 
small corrections in comparison to /i^. 

The requirement of cancellation of the small- a: singu- 
larities in the effective dynamics imposes some pertur- 
batively determined constraints on the otherwise non- 
perturbative ansatz for /i^. These constraints restrict 



the behavior of as a function of the gluon motion 
with respect to other constituents. In future, the refined 
versions of the same procedure may provide constraints 
that come closer to the actual behavior of gluons. This 
behavior is hoped to be uncovered in computer simu- 
lations that one may build around the first approxima- 
tion. Then, the extrapolation to ao = cus can recover the 
original theory from a few terms in the weak-coupling 
expansion if /z^ approximates the behavior of effective 
gluons well. Initially, the gluon mass term is viewed as a 
function of the relative momenta and Aqcd- The latter 
depends on ao as Aoexp(— c/ao) with a positive con- 
stant c. This means that Aqcd vanishes to all orders 
in the perturbative expansion, and /i^ is considered to 
be on the order of 1. Before one knows more, fi can 
only be estimated on the basis of implications for the 
resulting Schrodinger equation. The size of fi can be 
compared with four scales: Aqcd, Aq, m, and the Bohr 
momentum scale, fcs = aoTO/2, which is distinguished 
non-perturbatively by the Coulomb interaction. But if 
fi'^ is right, then [V — (ao/as)"/Lt^] must be a source of 
only small corrections in the whole range of couplings be- 
tween and as- This is taken for granted in the present 
article. Since the first approximation turns out to be 
not sensitive to the details of /i, new information can be 
obtained only in the refined calculations. 

The value of the weak coupling expansion scheme for 
Hamiltonians is that it starts from a local theory and 
leads to Hx that is capable of describing physically rele- 
vant non-perturbative dynamics even if Hx is calculated 
only in low orders. This idea is known to work in the case 
of QED: the Coulomb potential accounts for highly non- 
perturbative dynamics of atoms, including the nature of 
chemical bond, while the Hamiltonian itself is only of the 
formal order of a. Condensed matter physics illustrates 
this point in still wider domain. 

But when one applies the weak coupling expansion idea 
to QCD 2], one faces the fact that the strong coupling 
constant rises to 10, 30, or even 100 times larger value 
than in the case of atoms or positronium in QED. This 
leads to a complex interplay between the perturbative 
and non-perturbative parts of the calculation, enhanced 
dependence of observables on the RG parameter A, prob- 
lems with obtaining the Poincare symmetry in solutions, 
and amplification of artifacts due to the small-fc+ reg- 
ularization. Most of the problems seem to come from 
perturbation theory in the RG part of the calculation. 
An exact RG procedure by definition provides Hx whose 
structure depends on A but the spectra and 5'-matrix el- 
ements do not. However, when one uses expansion in 
powers of ao and then extrapolates to ao = a^ , a consid- 
erable dependence of the eigensolutions on A can ensue 
because of missing many terms. This is visible in mod- 
els that are asymptotically free and produce bound states 
jlEll^- One has to find the right value for as at given Ao 
from fits to bound-state observables, and perform consis- 
tency checks for whole sets of different observables SOj. 
In QCD, such checks involve the unknown functions of 



3 



momenta in the finite parts of the ultraviolet countert- 
erms, unknown terms depending on I^qcDi and artifacts 
of the regularization of small-Zc"*" divergences. So many 
unknowns suggest the possibility that the approach may 
never achieve the desired level of predictive power. But 
the simplicity of the harmonic potential found here in the 
first approximation illustrates that there is a high degree 
of order in the rich structure of H\. It is a consequence 
of preserving all seven kinematical symmetries of the LF 
scheme in the RG procedure. These symmetries limit the 
large number of terms that are allowed by the LF power 
counting using absolute momentum variables 0, to a 
much smaller number of terms that depend only on the 
relative momenta of the constituents. The LF symme- 
tries must also be respected by the initial regularization 
prescription for the required counterterms to be simple. 

The renormalization group procedure for effective par- 
ticles (RGEP), which is employed here, regulates the 
ultra-violet and small-fc'*' divergences by vertex factors 
r that are inserted only in the interaction terms. In the 
case of the small-fc"*" singularities, these factors limit only 
the ratios of momenta fc+. The ratios are limited from 
below by a dimensionless parameter 5. Details of the 
factors are explained later. Every creation or annihi- 
lation operator, labeled by momentum fc+, enters the ini- 
tial interaction Hamiltonian together with a correspond- 
ing factor rs(x), where x = /p~^ and p+ is the sum of 
all momenta that label all creation operators, or, equiva- 
lently, all annihilation operators in the same interaction 
term (see Appendix 0. In QCD with rs{x) ~ , H\ 
contains the coupling constant which depends on the 
scale A in the same asymptotically free way '3^ that 
characterizes the running coupling constant dependence 
on the renormalization scale in the Feynman diagrams 

Hamiltonians H\g with small Aq are worth studying be- 
cause their eigenstates can be expanded in the effective 
particle basis in the Fock space and the wave functions 
in this expansion are expected to correspond to the con- 
stituent picture of hadrons. This is envisioned in analogy 
to the models based on Yukawa theory 35]. The inter- 
actions are suppressed by the form factors and cannot 
copiously create new constituents, even if the coupling 
constant ao becomes large. Exotic hadron states may 
have their probability distributions shifted in the number 
of effective particles above the constituent quark model 
values of 2 or 3 [s^. The effective dynamics can be in 
agreement with requirements of special relativity even if 
it is limited to a small number of effective constituents, 
and the RGEP provides rules for constructing the rep- 
resentation of the Poincare group [s^- But the key fea- 
ture is that the transition between the bare and effective 
degrees of freedom is made in one and the same formal- 
ism. There is no need to match different formulations of 
the theory, such as in the case of lattice theory and the 
continuum perturbation theory in the Minkowski space 
HM 03 1 with none of the parts meant to cover the 
whole range of relevant scales on its own. 



One more comment is required concerning the small- 
X divergences in the eigenvalue equation for Hx- In the 
initial studies that used (5+ to limit bare particles' mo- 
menta and employed coupling coherence to derive certain 
H\ lile J2»^^i45i|, one could keep only a quark-anti- 
quark sector in the corresponding eigenvalue problem and 
the resulting equation was finite in the limit J"*" — > 0. 
Similarly, no small-Zc'^ divergences were encountered in 
the case of gluonium approximated by states of only two 
gluons In contrast, the present RGEP approach 

demands inclusion of states that contain an additional 
effective gluon which is needed to cancel small-x diver- 
gences. For example, if one keeps only a pair of the ef- 
fective quark and anti-quark, the leading non-relativistic 
(NR) terms are free from the small- a; singularities [s^ 
but relativistic corrections are singular j47, J£] . When 
the additional gluon is included, the condition of cancel- 
lation of the small- a; divergences becomes a guide in un- 
derstanding the gluon dynamics. The rules of including 
the gluons must be brought under quantitative control 
and the well-known case of heavy quarkonia provides a 
laboratory for testing the approach based on the gluon 
mass ansatz. The tests require a first approximation to 
begin with and a candidate is identified in the next sec- 
tions. 

This paper is organized as follows. SectionlTlldescribes 
the initial Hamiltonian of LF QGD with one heavy fla- 
vor and the procedure for deriving an effective Hx. Sec- 
tion mil discusses the eigenvalue equation for a quarko- 
nium, and introduces the ansatz for the gluon mass term. 
Small- a; effects in the dynamics are described in Section 
IIVI The resulting potential in the Schrodinger equation 
for a QQ bound state is described in Section Section 
IVII provides a brief summary and outlook. Appendices 
contain key details required for completeness. 



II. HAMILTONIANS 

The regularized canonical Hamiltonian of LF QCD 
with one heavy flavor of quarks, H , is given in Appendix 
IXI It includes ultraviolet counterterms. This section de- 
scribes the main features of H and the RGEP derivation 
of the effective Hamiltonian Hx with a finite width A. Hx 
is independent of the ultraviolet regularization factors 7'a 
in H when A ^ oo. The small- a; regularization factors 
rs, which are also present in H , and their role in Hx, will 
be discussed later. 

The initial Hamiltonian has the structure 

H = H^2 + Hj^2 + HjpAjf, + H(^^^Y + X (2) 

where the term H^i denotes the kinetic energy opera- 
tor for quarks, Ha2 the kinetic energy operator for glu- 
ons, H^Atp is the interaction term that couples gluons 
to quarks, i?(^^-)2 is the instantaneous interaction be- 
tween quarks, and X denotes all other terms including 
the counterterms. 
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In terms of the creation operators for bare particles, 
for quarks, for anti-quarks, for gluons, and the 
corresponding annihilation operators, the kinetic energy 
terms are of the form. 



and 



E 



k+ 



„bka 



Ha2 = 



ac 



'kl- 



,±2 



k+ 



^kcrc 



aha 



(3) 



(4) 



where k denotes the three kinematical momentum com- 
ponents, k~^ and k-^ = (k^, k^). The symbol in a bracket, 
such as [k], refers to the integration measure. 



[k] 



dk+d^k^ 
167r3fc+ 



(5) 



The subscript c stands for color and a for spin. The mass 
TO is assumed to be very large in comparison to Aqcd- 

The quark-gluon coupling terms in H.^Aip that preserve 
the number of quarks and anti-quarks, have the form 

F = g ^ / [123] 7=3,1 ^23 4«i^3 - j23 dla\d3 + h.c. 

123 

_ (6) 

The regularization factor 7=3,1 is singled out to indicate 
its presence. The coefficients J23 and j'23, are functions 
of the quark and giuon colors, spins, and momenta, with 
all details provided in the Appendix^ These coefficients 
contain the three-momentum conservation 5-function fac- 
tors, denoted by (5, color factors ^23: and products of 
spinors, ^23 = U2j'^U3 and ^'32 = W37'^W2- The latter are 
contracted with polarization vectors for gluons, so that 



J23 



^ ^23 dfJ-v ^23 ^1 



and 



J23 = S gp-v j32 ^1 ' 



(7) 



(8) 



The instantaneous term H^^^^i contains 
Z = -g^Y. I [1234]5t?2i23 [3t23tJ{kt 

T no A " 



+^2l 



kt) 



X ['^1,2^4,3 + ^2,1^3,4] b\dld4b2 



(9) 



The current factors j and the gluon polarization vec- 
tors grow with the relative transverse momenta of the 
interacting particles, k"*". These can increase to infinity 
and the regularization factors are introduced to limit 
the range to a finite one. In addition, there are small-x 
divergences due to the inverse powers of x, especially in 
the gluon polarization vectors that contain terms propor- 
tional to K'^/x. The small-x singularities are regulated 
by factors rg. 

The RGEP procedure generates ultraviolet countert- 
erms contained in the operator X in Eq. ^ and renders 



the effective particle Hamiltonian H\ which is indepen- 
dent of rA . The procedure is defined order by order in 
the formal expansion in powers of the bare coupling con- 
stant g. This expansion is eventually re- written in terms 
of the effective coupling constant, gx, which replaces g 
in H\ and depends on the ratio X/Aqcd 32]. The pro- 
cedure is designed so that energy differences in denom- 
inators of the perturbative evaluation of Hx are limited 
from below by A. Therefore, no infrared divergences are 
encountered in the evaluation of Hx. Also, no pertur- 
bative intrusion into the binding mechanism is generated 
when A is kept above the scale of typical relative momenta 
of the bound-state constituents. These features qualify 
the RGEP as a candidate for providing an answer to the 
well-known question of how it is possible that a simple 
two-body Schrodinger equation ma y re present a solution 
to a theory as complex as QCD |^l49j]. 

A very brief recapitulation of the RGEP is provided 
here for completeness. The derivation of Hx begins with 
a unitary change of the degrees of freedom from the bare 
quarks and gluons in Eq. ij^Jl to the effective ones. Let 
q commonly denote the bare operators b\ d^ , and a\ 
and their Hermitian conjugates. The operators q are 
transformed by a unitary operator Ux into operators qx 
that create or annihilate effective particles with identical 
quantum numbers. 



q\ 



Ux q Ul 



(10) 



The bare point- like particles in H of Eq. jSJ correspond 
to A equal infinity. One rewrites the Hamiltonian H in 
terms of qx and obtains 



H = Hx{qx) ■ 



Using Ux , one has 



Hx{q) 



uIhUx- 



(11) 



(12) 



Thus, Tix has the same coefficient functions in front of 
products of qs as the effective Hx has in front of the 
unitarily equivalent products of qx&. Differentiating Tix 
with respect to A, one obtains 



'H'x — — [Ta, Hx] 



(13) 



where Tx — Uj^W-^. 7^ is constructed using the notion of 
vertex form factors for effective particles. For example, 
if an operator without a form factor has the structure 



Ox= l[123]Vx{l,2,3)ql,ql^qx3 



(14) 



then the operator with a form factor is written as fxOx 
and has the structure 



fxOx = I [123] /a(Mi2,M3) 14(1,2,3) qWx2l>^3 



(15) 
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Different choices of the function f\ imply different inter- 
actions. The choice adopted in this study is 32J 



fxiMi2.M3) = eM~iMi,~Miy/X^)] 



(16) 



For any operator O expressible as a linear combination of 
products of creation and annihilation operators, fO con- 
tains a form factor f\{Mc, A^a) in front of every product. 
A4c and Ada stand for the total free invariant masses of 
the particles created (subscript c) and annihilated (sub- 
script a) by a given product. 

The effective Hamiltonian is defined to have the struc- 
ture 



^^A — f\G\ , 



(17) 



where G\ has to be calculated for given f\. One uses 
Q\ = G\{q), which is introduced in the same way as Tix 
in Eq. Ijlll) . Qi satisfies the differential equation 



Q'l = [fQi. {(1 - f)Ql}'g,] , 



(18) 



where Gi — G — G^), Go '^^ the part of H that does not 
depend on the coupling constant (/, and the curly bracket 
with subscript fJo denotes T that solves [sj 

[r,Go] = [{l-f)Gi]'. (19) 
The initial condition for Eq. (|18|) is that Goo — H, 

Gx = H + f ds [fsGls,{i'^'fs)Gls}g,] ■ (20) 
J oc 

This equation shows that one can find the counterterms 
X in H that remove regularization dependence from Gx ■ 
T^x = fxGx and Hx is obtained by replacing q hy qx- 
Gix is expanded into a series of terms t„ ^ g", 



Gi 



E 



(21) 



Ti is independent of A. Only the term H^a^ needs to be 
discussed here. According to Eq. ©, ti =a2i+cki2, 
where a2i denotes terms that create a gluon and ai2 the 
terms that annihilate a gluon (the left subscript denotes 
the number of creation and the right subscript the num- 
ber of annihilation operators). The corresponding effec- 
tive Hamiltonian interaction term is obtained by multi- 
plying the integrand in Eq. © by fx and transforming 
qs into qx^- 

When one neglects the terms that change the number 
of particles by more than one, T2 = /3ii -I- /322- Equation 
(|TS|l implies 



^2 = [{/'n},/Tl]=/2[TlTl] 



(22) 



with /2 = {/'}/-/{/'}. The first factor / in /2 refers to 
invariant masses in the first r in the square bracket, and 
the second / in /2 is for the second r. The square bracket 
denotes all connected terms that result from contractions 



that replace products gi^j by commutators [gi,(7j]. The 
solution for T2 is then given by 



T2X = ^2 A kin] 



T2c 



(23) 



where T200 includes 77(^^)2 and the second-order mass 
counterterms from X in Eq. |(2Jl. T2X depends on incom- 
ing and outgoing momenta in the two vertices generated 
by the ts. If one labels the three successive configura- 
tions of particle momenta by a, 6, and c, in the sequence 
aTabbTbcC, and introduces the symbol uv = A4^y — A^^j 
where A^J„ denotes the free invariant mass of a set of 
particles from the configuration u that are connected to 
the particles in the configuration v by the interaction Tuv 
in the sequence uTuvV, the vertex form factor of Eq. I|16|) 
in the interaction can be written as 

fx{Mab,Mba) = cxp[-(a6VA')] = fab- (24) 

If one then denotes the parent momentum for the vertex 
T„t, by Puvi and writes Puv in place of P^^, while all the 
minus components of momenta of the virtual quarks and 
gluons are given by the eigenvalues of Go = H^p^ + H^^ , 

, . Pbaba+pbcbc 
•^2 (a, b, c) = ^ [fabfbc - 1] • (25) 

The second-order perturbation theory renders 

Hx = TqX + TgX + fx [YqgX + V,gX + Zqqx] ■ (26) 

The kinetic energy term for effective quarks is 

,k^^ + nv 



where 



k+ 



Xkac 



Xkcrc 



dxkac 



(27) 



2 2 , 

mx = + 



(4/3) 5' f[xK]rUx)Y.\r2,e:, 
•' 12 



X [J^2A(™^A^^m2) - Jc■2Ao(™^^^^m2)] /^+_ 

(28) 

In the order of appearance, ttiq is the quark mass squared 
that should be present in Tqx^ in order to fit data for 
quarkonia, m§ — m? + o{g^). Factor 4/3 comes from 
color, {N"^ — 1)1 {2Nc)- The integration measure is. 



[xk] = dxd?K.^ /[IQn'^xil - x)] 



(29) 



where x — /k^ is the fraction of the quark momentum 
fca carried by the virtual gluon, and — kj- — xk^ is the 
relative transverse momentum of the gluon with respect 
to the quark 2. The effective mass does not depend on the 
particle motion. This is a unique property of the RGEP 
in LF dynamics. The small-x regularization factor is 



where 132 



fs{x) = rs{x) rs{l - x) , 



rsix) = x^0{x) 



(30) 



(31) 
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The middle argument of Ti\ is 

^{n? + n^'^)l{\-x) + K^'^lx, (32) 

and 

^2A(™^>^^m2)/fc3+ = {il{M\ m") - 1] l{M^ - ii?) . 

(33) 

The gluon kinetic energy term reads 



(34) 



The expUcit form of /i^ "sT, '35| is not needed here. 
The next term in Eq. (|^ is Y\ — f\Yqg\, 

= .9E123 /[123]r5(xi/3)r5(x2/3)/A(X?2,"^') 

X j23 ^L^Al^AS + j23 42«1Wa3 + ^-C . (35) 

The effective potential term, V\ — f\Vqq\, originates 
from the exchange of bare gluons with jumps in the in- 
variant mass of intermediate states above A. 



Vx 



where 



I [12M] 5 tl^tl^Vxili, 24.) h\dldih2. 

1234 



(36) 



^a(13,24) = 



X [0iz)fsix5/i)rs{x5/4) J^2a(1, 253, 4) 

(37) 




FIG. 1: Momentum labels in the interaction term mediated 
by the exchange of one high-energy gluon. The same labeling 
is used in the exchange of effective low-energy gluon in the 
next section and appendices. 

The sum over polarizations of the intermediate gluon 
reads 



(38) 



where the gluon momentum is 

k+^^=eiz){k+'^-k+'^), (39) 
and e(z) = e{z)-9{-z), 

z^{k+- k+)/{k+ + k+) , (40) 



while x^ = \z\ = k^ /{k^ + k'^), and 
^5 ~ ^t'^/^b ■ 



(41) 



The last term in Eq. ]2<i\i is the instantaneous interac- 
tion between effective quarks, Z\ = f\Zqq\. 



Zx = ~9^Y. f m4]K2tt3Zxil'i,24)bi,d{^dMbx2, 

1234 



(42) 



where 



1 



Za(13,24) = j + f^(Mi,,Mi^) 
X [d{z)fs{x5/i)fs{xr,/4) 

+ e{-z)Psix5/3)Ps{x5/2)]- (43) 



1^- 



3 > ' > 4 

FIG. 2: Momentum labels in the instantaneous gluon inter- 
action term. 



III. EIGENVALUE EQUATION 

Once A is lowered in perturbation theory to some value 
Ao just above the scale of binding mechanism, the re- 
sulting can produce the mass and wave function 
of a bound state of interest in a numerical diagonaliza- 
tion. The basis states can be limited to only those that 
have free invariant masses within a range of size about A 
around the eigenvalue. This has been verified numerically 
in a matrix model with asymptotic freedom and bound 
states |2^ |2^ |33| • In that calculation, H\ was derived 
using perturbation theory up to 6th order. A quite small 
set of effective basis states with energies between 4 MeV 
and 4 GeV was sufficient to reach accuracy close to 1% in 
the computation of the bound-state energy on the order 
of 1 GeV. In great contrast, the initial Hamiltonian of 
the model coupled all states in the entire range between 
0.5 KeV and 65 TeV. Preliminary estimates performed in 
Yukawa-like theories also indicate that the form factors 
fx suppress large momentum changes so strongly that 
the effective dynamics derived in low-order perturbation 
theory receives only small corrections from higher order 
terms, even when the coupling constant is made compa- 
rable with 1 js^. In the case of heavy quarkonia the 
same strategy should work even more accurately than 
in the Yukawa theory because as may be very small in 
comparison to 1. But the weak coupling expansion for 
Hxo produces new interaction terms already in order ao. 
These are derived here. An ansatz for the gluon mass 
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allows us then to finesse the structure of the first approx- 
imation for the resulting QQ-potential. 
The eigenvalue problem for H\ reads 

H^\P)^E\P), (44) 

where \P) denotes an eigenstate of the operators and 
P^ with their eigenvalues denoted by P^ and P^ (an 
example of the RGEP construction of the Poincare alge- 
bra at scale A in quantum field theory is given in [33 )■ 
The eigenvalue E has the form 

E = (M2 + P-L2)/p+. (45) 

The center-of-mass motion is separated from the bind- 
ing mechanism, which is a unique LF-dynamics feature 
preserved by the RGEP, and P~^ and P^ drop out of 
the eigenvalue equation. The state \P) is written in the 
effective particle basis as 

l^') = \QxQx) + \QxQx9x) + ■■■ ■ (46) 

The dots denote components with more than three ef- 
fective particles. Such expansion does not apply in the 
case of the bare particles because those interact locally 
and the interactions disperse probability densi ty t o high 
momentum regions and multi-particle sectors |35| . The 
wave function ipi3{K^^,xi) of the effective valence com- 
ponent \QxQx)i is introduced by the formula 

\QxQx) = J [rs]P+~S^Pu{Kis,xi)bl^dl^\o), (47) 

where the quark and anti-quark quantum numbers are la- 
beled with 1 and 3, respectively. ipi3{K^^,xi) must have 
dimension of l/^^j^g for the canonical normalization con- 
dition to give {P'\P) = P^S{P — P') (quantum numbers 
of the states \P) and \P') must be the same). The rel- 
ative transverse momentum of two particles, 1 and 3, is 
always defined as 

'^iz - {k+ki - k+ki)/ik+ + k+) , (48) 

and xi = k^ / P^ = 1— .T3. The wave function depends on 
A and quickly vanishes for > A. The normalization 
condition gives, 

{QxQx\QxQx) = Nqq{\) P+5{P - P')\p,=p , (49) 

where 

^QqW = E /[2^i'«i3]IV'i3(«^3,a;i)p. (50) 

13 •' 

The probability of finding other components than the 
\QxQx) is given by 1 — Nqq{\). The value of Nqq{(X)) is 
not known but it may be close to 0. On the other hand, 
one expects Nqq{X) to be close to 1 when m ^ Aqcd 
and 

to>A>Aqc_d. (51) 



When the wave function is negligible for relative mo- 
menta much larger than such A, the NR approximation 
must be accurate in description of the relative motion of 
quarks. In addition, when a\ 1, the Coulomb binding 
mechanism is expected to work and the dominant region 
of momenta should lie around the Bohr momentum scale 
ks — a\ m, provided that A 3> kp- At the same time, all 
the fermion spin and relativistic correction factors can- 
not become large (or diverging JsJ) in the NR expansion 
because of the presence of f\ [33 • But the dynamics of 
the dominant \Q\Q\) component receives some signifi- 
cant contributions from the \QxQxgx) component in the 
small- 0:5 region, since the coupling to the gluon sector 
grows like k^/x5 when x^ ^ 0. The gluon component 
may have a negligible contribution to the norm but has 
to be accounted for when 2:5 is small. 

If one neglected sectors with gluons entirely, the eigen- 
value Eq. 1)44(1 would read 

[Tqx + fx iV,,-x + Z,qx)] \QxQx) = E\QxQx) ■ (52) 

This equation is mentioned here because an analogous 
one was considered before [i^, |3 in a scheme us- 
ing coupling coherence and the absolute lower bound on 
gluon momenta, k'^ > . The equation found in "43] had 
a finite limit when (5+^0. The resulting dynamics con- 
tained a logarithmically rising potential and reproduced 
some of the characteristic features of the charmonium 
and bottomonium spectra. This was a considerable suc- 
cess in view of how crude were the approximations and 
the fact that the potential derived in and employed 
in 44, 45] behaved differently in the transverse and lon- 
gitudinal directions. But that strategy could not work in 
the RGEP approach. 

Three reasons can be given now for why the pure 
\QxQx) approximation is not allowed in solving the eigen- 
value problem for Hx. Two of them are related to the 
fact that the couphng between the sectors \Q\Q\) and 
\QxQxgx) is proportional to the first power of the cou- 
pling constant g, being mediated by the term Y\ of Eq. 
(|35|l . The first argument is non-perturbative. It is based 
on the result that the matrix models with asymptotic 
freedom and bound states lead to a successful approx- 
imation to the eigenvalue equation with H\ of second 
order in g\ if and only if all important matrix elements 
in the properly chosen energy window are accounted for. 
These certainly include matrix elements on the order of g 
[2^123. The second argument is perturbative, and con- 
cerns the evaluation of the effective Hamiltonian that acts 
in the sector \Q\Q\) alone. When the states \QxQxgx) 
are lifted in energy by an amount of order 1, quantum 
transitions in the quark-anti-quark sector that proceed 
through the states \QxQxgx) are formally of order g^ 
and must be included when one computes the quark-anti- 
quark dynamics in a series of powers of g up to terms of 
the explicit order of a\. The third argument is based on 
the fact that Eq. H52|) has a finite limit when 6^0 only 
in the leading NR approximatio n l3ll . Relativistic cor- 
rections contain singularities [47I |48| and the additional 
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gluon sector has to be taken into account to remove them. 

The eigenvalue Eq. imphes that the \QxQ\g\) 

component satisfies equation 



[Tqx + Tgx + ^msA - E] \QxQx9x) = -yx\QxQx 



(53) 



K/ijgA denotes interactions with sectors with more than 
one gluon and/or additional quark-anti-quark pairs, and 
the non-abelian gluon-quark and quark-anti-quark po- 
tentials of order ax- The interactions cause a shift in 
the gluon energy and make the eigenvalue equation differ 
from a similar one for positronium. The idea of seeking 
an ansatz for the shift and building a corresponding first 
approximation in the quark-anti-quark sector may seem 
completely new, but it patterns QED with the exception 
that there one knows from the outset that the leading ap- 
proximation to a Hydrogen atom or positronium is given 
by a two- body Schrodinger equation with a Coulomb po- 
tential [5^ IslL Is^ . The NR lattice approach to heavy 
quarkonia |53 . l54| also starts from a two-body picture. 
The key argument is not theoretical but comes from the 
phenomenology of hadrons. Theoretically, an ansatz for 
the energy shift in the sector \QxQxgx) is an attempt 
to harness the giant eigenvalue problem for H by turn- 
ing it into the eigenvalue problem for Hx and identifying 
the corrected Coulomb picture that may apply as a first 
approximation in QCD. 

A practical way to increase the invariant mass of the 
three-body sector and preserve the kinematical symme- 
tries of LF dynamics is to add a mass to Tgx, using 
the rules outlined in Section ^ see Eq. Since the 

divergence in the small- a; region disappears in the case of 
positronium when one adds a sector with a massless pho- 
ton, a gluon mass that approaches zero when 0:5 —> can 
remove the small-xs divergence in the quarkonium case. 
The rotational symmetry condition on /i^ can be imposed 
by demanding that the resulting potential in the quark- 
anti-quark sector is a rotationally symmetric function in 
the center-of-mass variables (see below) . 

Given a gluon mass ansatz, the whole eigenvalue prob- 
lem for Hx is limited to only two coupled equations. 



Y\QQg) 



fg)\QQg) + Y\QQ) = E\QQg) 
^f{V,^ + Z,^)]\QQ) = E\QQ) 



(54) 



(the subscript A is omitted). Tgx differs from Tgx of Eq. 
by replacement of the perturbative by the ansatz 
It is understood that /z^ may depend on A. All terms 
of order g^ in the three-body sector are ignored because 
they do not contribute to the dynamics of the \QQ) com- 
ponent in second-order perturbation theory (see below). 
This dynamics is described by the Hamiltonian Hqq that 
acts only in the quark-anti-quark sector. One should keep 
in Hqq terms of formal orders I, g, and g^, when the ef- 
fective Hamiltonian Hx is calculated up to the terms of 
order g'^, while /^^ ~ 1. The bare g is understood to 
go over in higher order calculations to a suitably defined 
gx 113, 112 • The perturbative expansion is applied only 
in the evaluation of Hqq. Solution for the bound-state 
spectrum of Hqq is not perturbative. 

To evaluate Hqq as a power series in g, one can intro- 
duce an operator R that expresses the 3-body component 
through the 2-body one. 



\QQg) = R\QQ) . 



(55) 



Since Y is of order g, R\s expected to be at least of order 
g. If V denotes the projector on the \QQ) sector, one has 
= (1 - p)i? = RV and 'PR_= R{1 - V) = 0. The 
effective Hamiltonian in the \QQ) sector is then given by 
the formula (see also ||35| concerning the context 

of RGEP) 



Hr 



1 



In the first order in g, 



{V + R'')Hx{V + R) 



VV + R^R 

(56) 



RTg-iTg+fg)R = Y. 



(57) 



Consequently, the second-order expression for the matrix 
elements of Hqq between different states i and j in the 
\QQ) sector, is 



{i\HQQ\j) = miTg +f{V,, + Z,,)] \j) 
( I 1 ^ 



Ej Tq Tg Ei Tq Tg 



y\3)- (58) 



The effective eigenvalue equation for heavy quarkonia, 
Hqq\QQ) = E\QQ), takes the form 



■-IS 



+ ml m^(l) m|,(3) 



X1X3 



Xl 



X3 



t(;{K^3,xi) - ^ 
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3 167r3 



X2Xi 



■vx{13,24)^{4^,X2) = 0, 



(59) 



where 



vxiU,24) = Fa(13,24) + Za(13,24) 



2X5 



d^uih) j'^^Jis w'a(13, 24) 



(60) 
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and 



i"a(13,24) 



+ m^)/x, - - ^2^-4)2 + ^2(2, 5, 3)1/2:5 



^24 



'■)/x2 



e{-z) fs{xz/3)fs{x5/2) .f\{m'^iMli)f\{Ml^,m^) 



24 



K^3)2 5, 4)]/x5 - (4^2 + ^2)/^^ 



6'(z) r5(x5/i)f5(x5/4) /a(to2^ 7V{2^)/^(7V{2_^^ ^2) 



+ m2)/a;4 - [(^5^3 - ^4)2 + /i2(2, 5, 3)]/a;5 - + m2)/x3 

^(-z) r^(a::5/3)fa(a;5/2) /A(m2, A^2^)/a(X^i, m2) 

(^^42 + m2)/x2 - [(^4 - K^3)2 + 5, 4)]/a;5 - (k^s' + m2)/a;i 



(61) 



The terms with 0(z) describe the emission of the gluon by the quark and absorption by the anti-quark, while the terms 
with 2:) describe the gluon emission by the anti-quark and absorption by the quark, see Fig. ^ The first square 
bracket corresponds to the first term in the large round bracket in Eq. (|58|l . and the second bracket corresponds to 
the second term. The mass terms, my, originate from the emission and re-absorption of the effective gluon by the 
same quark, in which case both terms in the bracket of Eq. (|58|l are equal. The mass terms read 



where 



m2,(l) = (4/3) g2 / {xK\r^,(^x) jl{m\M^) ■ 



{nf^ + rin?)/xi - {nf^ + MD/x 



(62) 



Ml = [n^^ + ^?{l',b',2,)]/x + {K^^ +m')l{l-x), 



(63) 



and, 



where 



m^(3) = (4/3) 5' / [xK]fi{x)fi{m\M') 



{Ki,^+m^)/x,~{4,^+Ml)/x, 



(64) 



Ml = [k^^ + ^i''{l,5',3')]/x + {k^^ + m'')/{l ^ x) . 



(65) 



M is given by Eq. H32() . The subscript 1' denotes the in- 
termediate quark and 5' denotes the intermediate gluon 
in the self- interaction of the effective quark 1, and, sim- 
ilarly, 3' and 5' denote the intermediate anti-quark and 
gluon in the self- interaction of the effective anti-quark 3. 

The gluon four-momentum fcs in the sum over polar- 
izations, i.e., in d^j^^{k^) in Eq. (|38|) . can be written as 

fcs" = £(z)g2. + [k^ ~ e{z)q7^ 12 , (66) 

where ij refers to quarks 1 and 2, or anti-quarks 3 and 4, 

(li] ^ h - kj . (67) 

The quark momentum four- vectors are on the mass shell. 
Since the gluon connects two vertices, one momentum k^ 
in d^j,,j{k^) is contracted with the current carried by the 
quark, and the other with the current of the anti-quark. 
In the self-interactions, both momenta are contracted 
with the same current. The momentum /cs contracted 
with current jfj can be expressed through qfy But the 
current conservation implies that the terms proportional 
to qij give zero. Therefore, one can replace Eq. (|38|) in 



the gluon exchange terms by 

dfivik^j = -g^ii^ + n^_,ny 
X [fcg- + e{z){k^ - fcf + fcg- - fc4-)/2] lk+ . (68) 

In the quark self-interaction one has 

dt^v{kv,,) = -g^^ + n^Ui, ^ f ^ — , (69) 

k+ 

with an analogous result for the anti-quark. 

The terms with the metric 5^,^ are regular in the small- 
X5 region, while the terms with n^^rii, are singular. The 
metric terms lead in the well-known way to the Breit- 
Fermi spin-dependent terms with a Coulomb potential. 
A discussion of the Breit-Fermi terms and gluons in the 
context of QCD can be found in ^tnd references 
therein. The singular terms with n^rii, are independent 
of the quark spin. It is shown below that the latter gen- 
erate the harmonic force between quarks when combined 
with the fermions' self-interactions, which are also in- 
dependent of the spin. Thus, the harmonic force ap- 
pears without Breit-Fermi terms. This result sets the 
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RGEP approach apart from the models employed in Refs. 
[ll|5l|5i|5lllfl.lllil|63, where one had to guess 
whether a confining potential appeared with or with- 
out Breit-Fermi terms. The spin-independent harmonic 
force is akin in this respect to the lattice picture and the 
orig inal charmonium model based on the Coulomb force 

[11 mil 113. 

In the explicit discussion of singular small- a; features of 
Eq. (|59() in the next section, all the g^i, terms are omit- 
ted. The reader should keep their presence in mind until 



they are re-inserted in Section The symbols of mass, 
wave function, and potential arc provided with a tilde 
as a reminder about the need to include the g^j^i, terms. 
Also, expressions for the quark masses are simplified by 
considering from now on only A = Aq. The subscript 
indicates that A = Aq. The ansatz for ^P' is understood 
to correspond to Aq. 

With the g^^ terms hidden and A = Aq , the eigenvalue 
equation reads 



/ nh? + ml •fni rh\ ~, 1 ^ 4 g^ /" dxodP'K.oA itoiti ^ , s ~, \ , 

" +^-f^-MM V(^j-3,Xi)---|-3 / 24 7^^ ~^(^3^24)^(44,X2)=0, (70) 

\ X1X3 x\ xz ) o loTT-' J X2X4 ^ xi^ 

where 

uo(13,24) = /i3,24 [^5" + e(2)(fc2" - fcf + ^3" - [^'(2)^5/17=5/4 J^i, 253,4 + ^(-^)^5/3^'5/2 -^^s, 154,2] 

+ /l3,24 [0{z)f;i/lf^/i + 6'(-z)f5/3f5^2] + ^ [^^ + £{z){k2 - + - V)/2] Wq , (71) 

and the last factor of wa = P^w\g{13, 24) is abbreviated to 

^(2)?'5/1^5/4/l,52/53,4 , ^'(-2) ?5/3^5/2 /3,54/51,2 , ^'(2) 7=5/17=5/4 71^52/53,4 z) f 5/3^5/2 13,54^1,2 



Wo 



k\ kg (2,5,3) ^2 k^ fcg (1,5,4) fcg (2,5,3) k^^ k2 fcg (1,5,4) k-^ 



(72) 



The compact notation includes 

U,,^h,XMlM]), (73) 
7=5/,, = fs{xz/i) , (74) 
^i^kj ^ J^2Xo{hk,j) , (75) 
fc5-(7, J, k) - [4^ + A^2(7, J, k)]/k+ , (76) 
4 = £iz) («^3 - 4i) ■ (77) 
The mass terms with the g^i, terms suppressed are 
~2 4 2 /• 2 2 + 1 x{M'-m^) 

(78) 

for z = 1,3. Ml is given by Eq. and M3 by Eq. 

(|()5|l . In both cases, is given by Eq. and the 

factor {j+lk+f = 4(1 - x). 

IV. SMALL- X BEHAVIOR 

All small- a; singularities of the eigenvalue Eq. (|59|l are 
contained in Eqs. lO, ((TIJ, and lO. We first discuss 
the exchange terms, then the mass terms, and finally the 
net effect of the interplay between these terms. 

The analysis hinges on the properties of the energy 
of motion of a gluon with respect to the parent quark. 



Pi" krj — ^1^5^ 

jx^. Pi is the momentum of the parent 
quark i. The momentum ranges under the integrals 
from to cx), while 3:5 can reach (in the mass terms, the 
integrals are over and x). Appendices IbI and O pro- 
vide definitions of all variables used in the description of 
the integrands. The key difficulty is that the ratio of two 
variables of different kinds, and x, varies quickly with 
a change of any one of them. This complexity is related 
to the power counting rules for the Hamiltonian densities 
on the LF But the analysis described here concerns 
only the relative motion of the effective particles and it is 
simplified by taking advantage of the NR limit after the 
finitcness of the small- a; dynamics with the gluon mass 
ansatz is established. 

The singularity in the effective gluon exchange term 
is tempered by the product of two vertex form factors, 
//. The form factors vanish exponentially fast when 
K^'^/x^ — > 00. This prevents 3:5 from becoming small 
unless vanishes at least as fast as Therefore, 
the measure of integration over transverse momenta is 
on the order of X5 when X5 and it reduces the di- 
vergence to a logarithmic one. The logarithmic diver- 
gence is taken care of using the gluon mass ansatz. The 
mechanism of reducing singularities to only logarithmic 
ones does not work in the instantaneous interaction term 
^^(13, 24) and in the terms in Vx(13, 24) that come with- 
out // in J^2X- But all the terms without // are inde- 
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pendent of fi^ and the dx^/x\ and dx^/x^ singularities 
cancel out in them perturbatively |l6j| . 

In the fermion self-interactions an analogous pattern 
of the singularities occurs. But one has to also consider 
the size of mp. The latter is determined by the size of 
the free ultraviolet-finite part of the quark mass coun- 
terterm in the initial Hamiltonian of Eq. ([SJl . That size 
is related to an ansatz for a gluon mass term in the sec- 
tors \Qg) and \Qg) in the eigenvalue equations for states 
with quantum numbers of a single fermion, see Appendix 
lO Eventually, the gluon mass ansatz leads to the re- 
sult that the single-quark eigenvalue diverges logarithmi- 
cally in the limit i5 — > 0, while the quark self-interaction 
in the quarkonium dynamics becomes finite. The self- 
interaction and effective gluon exchange, both finite due 
to the chosen behavior of the gluon mass ansatz, lead 
together to the harmonic potential which is described in 
the next section. 

According to the Appendix ^ the dominant exchange 
terms in Eq. (|70l) read 



{io(13,24) = B{z)i+iow + e{-z)i^io 



where v+iow is given in Eq. 
(|B20|) . In the hmit x^ 0, 



mm. and V- 



low 



— /l,52/l 



1,52^/53,4 ^2 



^2(2,5,3) 



^2(2,5,3) 



and 



V-lo 



73,54/51,2 



m'(1,5,4) 



(79) 
in Eq. 



(80) 



(81) 



Since is on the order of \z\, one obtains the result 
that if fj,^ vanishes faster than g^^, i.e., faster than X5, 
the potential produces a finite effect in the limit of (5 — ^ 0. 
In the denominators of Eqs. (|80() and (|81|l there also ap- 
pears — {2mz)^, which is negligible in comparison to 
the leading terms on the order of z but can be included 
here on the basis of hindsight to take advantage of the 
NR nature of the quarks' motion with respect to each 
other. The larger is the quark mass m for fixed Aq and 
the smaller is Aqcd, the more accurate the NR picture 
actually becomes after the small- a; divergences are re- 
moved. Writing = qt, with q = \q\, t = cos9, the 
singular factor l/xf equals Am^/ql = {^m? / q^)t~'^ . The 
integration measure d'^q is proportional to q^ and the 
small-x singularity is actually produced by the angle in- 
tegration dt/t^. fi^ should vanish for i ^ in order to 
remove the singularity. An example of such behavior is 
used below to provide a constructive context for the steps 
that follow. The final result is not sensitive to the details 
of the example. Given that /i^ vanishes faster than q^, 
one can write 



M^(«,5,j) = c^{i,5,j)q^ , 



(82) 



and determine behavior of c(i,5,j) from the condition 
that 



c(i,5, j) 



c^(»,5,j) 
l + c2(i,5,j) 



(83) 



should vanish for x^ ^ 0. 

The only information about the three-particle sector 
that is available in the relativistic construction of q and 
c(«, 5, j) are the _L and -I- components of the momenta ki, 

, and kj . Two physical constraints are used in defining 
a helpful vector q: the definition must respect all kine- 
matical LF symmetries (to render a boost-invariant de- 
scription of quarkonia), and it must reduce to qz = 2mz 

for z — !■ when |fci3|/m and |fc24|/m approach in the 
NR limit. A geometrically motivated candidate for q 
is provided by the difference between the square of the 
free invariant mass of three effective particles in the state 
I QQg) and the square of the invariant mass of the QQ- 
pair in this state. The difference reads 



,-L2 



Ml, 



M ■ = 



2:5(1 - 2^5) 



(84) 



Multiplication by a;5(l — X5) produces an expression that 
tends in the limit of X5 — > to the three-momentum 
transfer squared that appears in the energy denomina- 
tors in the small-x dynamics. The components of q are 
therefore defined as q-^ = and 



qz = zM 



I] ■ 



(85) 



Further analysis of all exchange terms shows that if the 
ansatz mass fj?' behaves like 



(86) 



low 



of 



(or like q'^x^'^) with < < 1/2, the factors v± 

Eqs. HIjl9|l and (|lj20|l vanish in the limit X5 — > as x^" 

independently of the terms in the energy denominators 
3/2 

on the order of or smaller. Thus, the gluon exchange 
term becomes finite when 



(87) 



c(j,5,j) = c{t) 

and c(t) is a function that behaves as 

c{t) ^ c\t\'^/\ 

for t — !• 0, with c a constant. 

With this ansatz, the quark mass terms also become 
finite in the limit 6^0. Appendix O shows details of 
how ttiq is chosen in agreement with the physical picture 
explained in the Introduction and at the beginning of this 
section. Eq. (|78|l gives 111^+ rhf = m? + 5mf with i = 1, 3 



and 



4ao 



<5m^ = IP j dxd'n^fl{x)fl{^\M') 
M'^-m'^ f 1 1 



TWq — -m? Mj — nn? 

(89) 

The function A^q is given by Eq. (jC9p in terms of the 
gluon mass function /Xq that must satisfy the condition 
noil . The quark self-energies are positive if 



(90) 
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A simple way to satisfy this condition is to set /ig = 0. 
Then, 



IJ.'^{i\5',j) + (k-Lz +x2m2)/(l - a;) 



(91) 



where i = 1 and j = 3 or vice versa. The factors rg are 
no longer needed. 

The small- a; regularization disappears from the quarko- 
nium dynamics entirely. Finite phenomenological param- 
eters that describe the small- a; behavior of the gluon mass 
ansatz, such as (5^ in Eq. (|86(l . become responsible for the 
regularization of the exchange and self-interaction terms, 
preserving their distinct properties. The mass terms grow 
when decreases, while the effective gluon exchange po- 
tential provides a negative contribution that increases in 
magnitude at similar rate and compensates the size of 
the masses at small momentum transfers. The net result 
is described in the next section. 



V. THE QQ SCHRODINGER EQUATION 

The condition (|51|) validates the NR and weak cou- 
pling limits after the small- a; divergences are removed by 
the gluon mass ansatz. This section then identifies the 
leading structure in Hqq in formal order of a^. The ad- 
ditional simplification in the case of small is that the 
dominant interaction in Eq. H59() becomes equal to the 
well-known Coulomb term and one can find the leading 
correction analytically. 

Equation (|59|) can be re-written using the relative 
three- momentum variables described in Appendix IbI see 
Eq. (|lj2|l . The integration measure is 



X2X4 M24 
and Eq. (|59|l takes the form 



Xl X3 



(92) 



^'(fcis) 



(27r)V"i^ + fc| 



:C/(fcl3,fc24)^(fc24) =0. 



24 



(93) 



The mass corrections include now the g^^^, terms that were 
suppressed in the previous section, (5m| = Smf + Srn^, 
and the potential is 



Ay/XiX2X3X4: s . , . . i 



where Vg denotes the g^i, contribution in the exchange 
term. 

Since the form factor /i3,24 cuts off changes of the rela- 
tive momenta above Aq exponentially fast, one can focus 
on the eigenstates with lowest IVP and take advantage 
of the conditions {ki^l <ti m and |/c24| ^ that are 
satisfied in the entire domain of physically relevant prob- 
ability distribution. For such states, one can expand Eq. 
(|93|l in powers of fc/m, with the exception of the form 
factors that are needed for convergence. The Coulomb 
force defines the momentum scale of the inverse of the 
quarkonium Bohr radius, fcs = r^^ = aQm/2. When 
Ao ^ kB, the form factor f 13^24 does not differ from 1 
in the dynamically dominant region; /i3,24 matters only 
when one extends the expansion to high powers of k/m. 
These would lead to divergent integrals with Coulomb 
wave functions and require counterterms 50, 51] . The 
latter are not needed here and the lowest terms domi- 
nate I35I. The binding energy, B, is small in comparison 
to m. Writing the quarkonium mass as M = 2m + B and 
neglecting ~ B^/m, one obtains 



TO 



B 



Smf 



2m 



6ml 
2m 



i>{k- 



13^ 



-I- 



(27r)3 



A:24) ^"(^24) = . 



(95) 



The structure of Vqq and the size of the mass corrections 
5m\ and 5to§ need explanation. 

The two vertex form factors that appear inside the 
exchange and mass terms in Eq. H59|) . have arguments 
given in Appendix IbI in Eqs. (jB8|) to HB13p . When one 
writes the product of the two vertex form factors in the 



NR limit as exp (- 



reads 



= V2™i^. 
Aq t Ao 



(96) 



The limit to/Aq ^ 1 enforces g ^ Aq, the more so the 
smaller is t. The Coulomb binding mechanism is intact 
for Ao as small as several times ks [35j |. which is much 
smaller than to in the weak coupling limit. Thus, the mo- 
mentum transfer q is much smaller than ks in all terms 
that contain //. These terms become then negligible in 
comparison to the Coulomb term, unless they have sin- 
gularly small denominator factors for small t. That is the 
case for the mass and exchange terms when 5^ becomes 
small. In the presence of the g^i, contributions that were 
omitted in Section Hvl these terms are found as follows. 
The Hamiltonian Hqq has the structure 

Hqq = + Sm'^ if 1,9 + ^,0) -6m'^ iff, g + n,fi) 



+ f{l- ff)[{g + n,0) + z] 
+ f{ff)[{g + n,fi) + z], 



(97) 



(94) 



where g denotes the g^i, terms, n denotes the singular 
Ufj^ni, terms, and z denotes the instantaneous terms. The 
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gluon mass ansatz in energy denominators is indicated 
by an extra variable in the brackets, and says that the 
gluon mass is 0. The last two terms can be re-arranged 
as 



W is large and negative when 5^ is small. The ex- 
change term tends to compensate the positive contribu- 
tion of the mass terms. This can be made transparent 
by re-writing Eq. H95|) as 



f(l-ff)[{g,0) + {n,0) + z] + 



(98) 



The contribution of (n, 0) + z in the first term vanishes in 
the leading NR limit, see Appendix IbI Two of the terms 
with {g, 0) combine to f{g, 0), and reduce to the Coulomb 
term with the Breit-Fermi spin corrections. The remain- 
ing terms, with / in front also being equivalent to 1, 



/(//)[(5,/i)-(g,0) + (n,/i)+z] 



(99) 



add to the Coulomb term and produce together Vqq in 
Eq. (|101|l . The mass terms can be re- written, in the 
same fashion, as 

if.f)Sm^[{g, 0) - (5, /i) + (n, 0) - {n, . (100) 

Expressions H99(l and (|100|l show that the exchange po- 
tential and the mass terms have similar structures with 
opposite signs. A change of variables from x and to 
Qz = xm and ~ in the mass terms produces inte- 
grals in which the factor // ensures that q = \q \ m 
and one can again use the expansion in powers of the 
ratio of q/m. Since the integrands are symmetric func- 
tions of qz, one can extend the integration to negative 
qz and divide the result by 2, which produces the same 
integrands as in the exchange terms. Hence, 

VQQ{q) - {l + BF)Vc{q) +W{q), (101) 

where BF denotes the Breit-Fermi spin-dependent fac- 
tors. 



4 47ra 



W{q) 



3 

X exp 



■ Ana 



mq 



1 

2 



+ q' 



(102) 



(103) 



with = 6{z) fi^{2, 5, 3) + e{-z) 5, 4), and 
(5m? f (Pc 



^q 



(27r)= 



W{q) , 



(104) 



with ^2 equal ^ji'^{V ,i) for i = 1 and ^I'^il.b' for 
i = 3. 

If the gluon mass ansatz is 0, = and the quarko- 
nium dynamics reduces to the same as in QED with addi- 
tional color charge factor 4/3. A finite gluon mass ansatz 
introduces new dynamics which is discussed in the re- 
maining part of this section. 



fc2 
TO 



B 



Sq 
(2^ 



{l + BF)Vc{q)iAk~q) 



(Pq 
(2^ 



W{q) tp{k 



(105) 



There is no need to trace the small relativistic correc- 
tions before the main NR picture is identified. Only this 
picture is discussed below. 

Since \q \ m W is constrained to values much smaller 
than fcs, one can expand the wave function in the 
Coulomb region under the integral in the Taylor series 
and consider the lowest terms as candidates for the first 
approximation, 



ij}{k - q) 



d ^ 
oh 



2 dh dk, 



(106) 



The terms with odd powers of q average to 0. The bilinear 
terms contain q^ times (1 — t^) times cos^ (j), or sin^ <j), for 
i — j — 1, 2, respectively, and t^, for i 
over (j) produces tt times a vector 



3. The integral 



w{t) = {l~t^,l~t^,2t^). 



The variable q can be changed to u of Eq. H96|) 
introducing the constant 

V2m 



^0 



one obtains the vector 



dtt{l-t^)w{t)T{t) , 



du — — — -r e 



that appears in the resulting interaction term: 



W, 



QQ 



3 2tt 



3 



dkf 



(107) 
, and 

(108) 

(109) 
(110) 

(111) 



The next non-vanishing terms in the Taylor expansion 
contain the fourth and higher even powers of q. They are 
expected to be small in the momentum region dominated 
by the Coulomb dynamics and do not count around the 
bottom of the harmonic potential. The remaining ques- 
tion is if the harmonic approximation can be rotationally 
symmetric. 
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The interaction Wqq given by Eq. (|lll|l is rotationally 
symmetric when all components of r are equal, or 



dtt{l-t^)il-'St^)T{t) = 0. 



(112) 



The function T(t) depends on /i in a limited way because 
the integral over u in Eq. (|ll()|l extends only from to 
about 1, b/t is large, and bfi/t » 1 produces T{t) — (3, 



/3 = V^/4. 



(113) 



Behavior of r(t) near t = does not matter because of 
the factor t in Eq. H112|) . and the condition (|86|l is of 
little consequence if raises quickly from at t = 0. 
For any ansatz of Eq. 188|l with a small 5^, 



T{t) 



cm 



f3, 



(114) 



which is equivalent to the constant (3 if c{t) ^ 1 for t 7^ 0. 
The factorization feature is independent of the shape of 
the RGEP form factor / and provides an opportunity to 
fit analytically to satisfy Eq. (|112|) . Suppose that Eq. 
is valid and 



cHt) 



c-" A (1 - pr) . 



l + c^t) 
Eq. (|112|l is satisfied when 

and 

t^^' (1 - pt^) 



1 _ c2 (1 - pt^) 



(115) 



(116) 



(117) 



leads to a rotationally symmetric harmonic oscillator po- 
tential. All components of r are equal f , 



f = (3c 
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(118) 



p varies between and 17/25 when (5^ varies between 
and 1/2 in accord with Eq. (|86|l . In the limit 
0, p ^ and c^it) ^ c^/{l - c^) outside the area of 
very small t. In this example, the size of the coefficient 
functions c(i,5,j) depends on the size of the constant c 
and grows to large values when c — > 1 . 

Every gluon mass ansatz, that is large outside the 
area of small t and vanishes abruptly for t — > 0, leads 
to a spherically symmetric harmonic oscillator potential 
with the spring constant 



4 a 

3 TT 



f, 



and 



f - /3/6. 



(119) 



(120) 



The number /3 depends on the shape of the RGEP form 
factor /, which reflects the dependence of the effective 
dynamics on the RG scheme, but (3 is stable for fabS of 
similar shapes as functions of ab/X^, see Eq. (|24|) . 

The first approximation for heavy quarkonium dynam- 
ics in position space can be defined by the Fourier trans- 
form of the eigenvalue equation for Hqq, with 



(r |fc ) = exp (ikf) . 



(121) 



This transform exists only in the relative motion vari- 
ables, since the motion of the quarkonium as a whole 
is described in a relativistic fashion and the relation be- 
tween the relative motion of quarks and the motion of the 
bound state as a whole does not coincide for large speeds 
with the one known in NR theory. The Schrodinger equa- 
tion reads 



2™ 



A, 



4a 

T 



BF 



k 



where k = mcj^/2, and 



4 a 
3 n 



X 



6%/2 



(122) 
(123) 



The number [(3/{6V2)]^/^ gives (7r/1152)i/4 - 0.23, 
which is large enough for the frequency uj to fit into the 
ball park of phenomenologically plausible scales when one 
allows sufficiently large A and a with some choice for m. 
The observed spectrum of charmonium is known to have 
an intermediate character between the Coulomb and har- 
monic oscillator spectra. But the key problem is to de- 
termine the size and direction of corrections that need 
to be included in order to compare the theory with data. 
The inclusion of light quarks requires quantitative under- 
standing of the mechanism of chiral symmetry breaking 
in the effective particle approach and a comprehensive 
phenomenological analysis demands further advance in 
the theory. 



VI. CONCLUSION 

The Coulomb interaction between quarks in heavy 
quarkonia is corrected by the potential well that is ex- 
cavated by the one effective gluon exchange in the over- 
lapping self-interaction gluon clouds of the quarks. At 
the bottom, the well shape is a quadratic function of the 
distance between the quarks. The resulting harmonic os- 
cillator force plays the role of a confining one in a limited 
range. At distances much larger than the Bohr radius 
the quadratic approximation stops working. Emission 
of additional gluons and pairs of light quarks will fur- 
ther change the rate of growth of the potential. The size 
of these effects should be computable in the present ap- 
proach by evaluating effective Hamiltonians order by or- 
der in a weak coupling expansion and solving eigenvalue 
problems for them numerically. 
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The effective particle approach is of interest because it 
describes the relative motion of quarks independently of 
the speed of the quarkonium as a whole. This result is ob- 
tained at the price of setting up QCD in its Hamiltonian 
version in LF dynamics, with a host of difhculties in the 
renormalization program that had to be overcome. Fur- 
ther advances in the RGEP and methods of solving the 
eigenvalue equations for Hamiltonians H\ are expected 
to reflect the well known features of interactions of rel- 
ativistic particles. The first approximation for Hqq can 
be expected to work well in the refined calculations be- 
cause it appears to be largely independent in its structure 
from the details of the RGEP vertex form factors and the 
gluon mass ansatz. The first approximation also appears 
to involve the least possible degree of complexity as a ba- 
sis around which a meaningful successive approximation 
scheme can emerge. A few percent accuracy in evalu- 
ating effective Hamiltonians is known to be achievable 
using essentially the same method in the case of elemen- 
tary matrix models with asymptotic freedom and bound 
states. 

Since the approach developed here is boost invariant, 
it can connect physical images of hadrons in different 
frames as soon as the hadron dynamics is understood in 
one of them. Although the light quarks are expected to 
behave differently from the heavy ones, one should note 
that the Schrodinger equation with H\ does not lead to 
the spread of probability towards large relative momenta 
and large numbers of effective particles. The spread is 
halted because the interaction terms in H\ contain form 
factors. These form factors are the reason for hope that 
the effective particle expansion may converge. 

Aside from QCD, the same scheme for setting up and 
solving quantum field theory should be tested in the case 
of QED. There, the effective mass ansatz for virtual pho- 
tons is much more restricted and small- a; effects are of 
less significance. On the other hand, QED is not asymp- 
totically free and its effective nature requires better un- 
derstanding. The RGEP approach may help in defining 
QED as an effective theory. But one needs to first ver- 
ify if perturbation theory with H\ can produce covariant 
S'-matrix in QED in orders higher than second. 



APPENDIX A: REGULARIZED LF 
HAMILTONIAN OF QCD 

The canonical LF Hamiltonians of gauge theories, sim- 
ilar to the Hamiltonians in the infinite momentum frame 
[68, ,69, 70] , are well known JX, , and extensive liter- 
ature exists on the light-like axial gauges [Zilzi. The 
Hamiltonian given below is further specified by inclusion 
of the ultraviolet and small- a; regularization factors that 
render a computable operator. This means that H does 
not require a separate regularization prescription for eval- 
uating loop i nteg rals. The same regularized Hamiltonian 
was used in jS^I but the quark terms needed here were 
not explicitly given there. H is supplied with countert- 



erms H/^g. Their structure is known from considerations 
similar to ^ . Details can be calculated using RGEP. The 
initial Lagrangian is 



with one flavor of quarks of mass m. 



4 tJ-v ' 



and 



D,. 



(Al) 



(A2) 



(A3) 



where A = A'^t", [t'',t''] = and Tr{tH^) = 5"^ /2. 

The classical Nether generator of evolution in a:+ takes 
the form (the Gauss law constraint is formally solved in 
yl+ = gauge and the counterterms are added as the last 
term from hindsight), 

H = H^2 + + Hji^a + H^i + HjpAjp 

+ Has, (A4) 
where each term is an integral over the LF hyperplane, 



and. 



n 



T~l-ipAip 

'H^AAap 

T^[dAA\'2 

[aAA](vv) 



Hi = / dx^cfx^Hi 



V [A^,ApnA'-,A^]\ 



^\AHd^fA\ 

gid^Al[A",A^]\ 
1 

\g^ [^^+A^,A^fJr^[^^+A^,A^]\ 

(A12) 
(A13) 



(A5) 

(A6) 

(A7) 
(A8) 
(A9) 
(AlO) 
(All) 



n 



1 



1 



:[id+ A-^ , A^]" 



(A14) 



A quantum Hamiltonian is introduced by expanding the 
fields into Fourier components at = and imposing 
commutation relations on the latter. They define cre- 
ation and annihilation operators for bare particles. 



7 —ikx I jt ikx 



(A15) 
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The integration measure is 

9{k+)k+ <P 



[k] 



167r3fc+ 



(A16) 



= 16TT^k+S^'Jcc'S^ik-k'). (A17) 

d^{k - k') = 6{k+ ~ k'+)S{k^ - k'^)S{k^ - k'^). The cre- 
ation and annihilation operators have the power-counting 
dimension 1/fc^ (the same result holds for gluons, see 
below). The spinors are given by Uku = B(k,m)ua-, and 
Vka — B{k,m)va-, where = Cu* = i7^u*. and 
are the spinors for the fermions at rest in the arbitrarily 
chosen frame of reference where the quantization proce- 
dure is introduced. The matrix B{k, m) represents the 
LF boost that turns a particle with mass m at rest to a 
moving one that has the momentum k, fc^ — rri^, 



B(k,m) 



1 



[A+/c++A_(TO + fc^a-^)], (A18) 



where A± = 707 /2. This matrix mixes fc+ with m and 
k^ . But the second term, of the type ^Jk-^ /k+ when one 
counts m and k^ as similar, results only from writing the 
interaction terms in a way that is short and convenient in 
calculations. The independent degrees of freedom, = 
A-i^Vj contain only the parts proportional to ^k^ jm. 
Thus, ?/'+ has the dimension of which in the 

position space on the LF leads to {x^\fxr^ ^] \ih ^ 
1/k^. The spinors at rest are. 



y2m 
Va = \/2m 







(A19) 
(A20) 



where — — 1(72X0- — fX-o-- The gluon field at x+ — 



and the commutation relations read. 



(A21) 



= 16TT^k+S„,Jcc,5^{k - k') . (A22) 



a and have the dimension 1/k-^. The polarization four- 
vectors are introduced using LF boosts as for fermions, 

- (4. = 0,e^, = 2fc^4/fc+,e^) , (A23) 

except that the boosts are applied to the polarization 
vectors ej^ — (0,0, e^) that correspond to the selected 
state of a gluon moving along the z-axis |l2j . Terms that 



contain the ratio k-^/k~^, which mix the transverse and 
longitudinal momenta, are again only a shorthand nota- 
tion for writing interactions. The independent transverse 
field A-^ contains only polarization vectors that have 



dimension 1. A has dimension of k 



which matches 



the required on the LF Q when a ^ 1/fc^, as 

promised. 

The kinetic energy operator for quarks, i?^2, is given 
in Eq. (j2Jl , and for gluons, iJ^2 , is given in Eq. Q . The 
triple-gluon interaction reads 

Ha^ - /[123]%^-p)fA5(3,l) 

123 

X ^gYi23a\ala3 + gY*23ala2ai . (A24) 

The symbols introduced in this operator occur in all other 
terms and require explanation for completeness, see also 
p}2.] . The conservation of momentum in the interaction 
vertices is enforced by the factor 



The regularization factors are given by 
where 



(A25) 



(A26) 



(A27) 



The symbol p refers to the parent momentum (half of the 
sum of all momenta of all particles coupled in a vertex) , 
and d to the daughter particle momentum, i.e., the mo- 
mentum of the particle emitted or absorbed in the vertex. 
The arguments of the regularization factors are defined 
by 



•^d/p^p 



^d/p — ^d /^p — ■^d/Xp . 

The functions used here are Is^l 

rA(«^') - exph^^VA^], 
rs{x) ~ 6{x — e)x^ , 



(A28) 
(A29) 



(A30) 
(A31) 



and e/S tends to when (5 — > 0. The gluon spin vertex 
factor reads 



123 



(A32) 



X2/3 



6263 ■ e^K- 



The simplified notation means: e = 
quartic gluon vertex is 



a^i/3 
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Ha'^ = ^ / [1234]5(p^ ~ T i234ai44'^4 + ^A* i2Siaia\a3a4, + i2344«3a2ai , (A33) 

1234"' 



where 



SA4 1234 = -[h+2,ifi,i{e\el-elei-e\ei-elel)r'^-r'^* 

+ h+3,ih.2 (ele*^ ■ £3^4 - £^£4 • slel) f--^-^ f----^ 
+ r3+2,3r4,i {elel ■ e^e^ - £*£4 • e^ej) jacic^j ^ 

^A* 1234 = ri+2,lf3+4,3 (£^£3 ' £^£4 - £^£4 ' £^£3) r^^^ T^^^ 

+ [^34^=2,4 + ri,3?4,2](£t£2 ' £3£4 - £^£4 • £2£3) r^^T--' 

+ [^3,2^1,4 + ?2,3?4,l] (£^£2 ' £^^4 " £^£3 ' £2£4) ^^-^ ^^'^^ 



(A34) 



(A35) 



The abbreviated notation for the regularization factors is fp^a = fAsiP^d)- Denoting t^j = xL^^Xjo the quark-gluon 
coupling is given by 



123 ■' 

where the spin vertex factors are 

U2 flUs. = \/x:i/x2 xl 

f\02 = V X3/X2 



U2 ^1^3 ^23 bla\b3 — V3 ^\V2 ^32 djaids + Ui ^3V2 &W2"3 + h.c. 



3 ;,tjt. 



X\ X3 



V3 ?-1^2 — V -^3/^2 <,_3 
Ul ^3V2 = \/x3/xi\/x3/x2 Xl 

The instantaneous fermion interaction reads 



-l(^Kjy3 X £3 ; H Ki/3£3 cr - TOiCr £3 

"^3 



X3 , 
C-2: 
€-2. 



^^VAAV; = / [1234]5(p+ - p) (9^2) • {} 

1234 



(A36) 

(A37) 
(A38) 
(A39) 

(A40) 



where the curly brackets {} contain the operators ordered according to the rule adb; 



{} = n+2,1^3+4,3 

+ [^,2^4,3 + r'2, 1^3,4] 
^3,4^1+2,1 til 



til ^^^^ b\ala3b, + til ^^ii^^44«2di 



Xl + X2 



Xl + X2 



b\ala2b, + til ^4^^^44a3rfi 



Xl — X2 



Xl — X2 



+ 



+ 



+ 



23 Xlr3](73e_4 ,t.t„t 



Xl + X2 



61^402013 + h.c. 



f2,ir3+4,3 tf, Xim'y^^^i-. ^t^t^t + 



2 y [n,2r3,4 + r2,ir4,3] t 

1 
2 

1 



23 



xl[23]a3g-4 

Xl — X2 



+ [n,3r2,4 + r3,ir4,2] i?4 



Xl — X2 

xl[32](73g_ 

Xl - X3 



1 ;>twt 



61^40203 + h.c. 



^ ~ .23 xt[2*3*]x4 , - - .32 XI[3*2*]X4 ^ ,t t t. , ^ 

r-4,3r-n-2,i , h ^4,2^+3,1 ql , — 6|a;^a;^64 + h.c. 



7: I n,2r3+4,4i?4 



Xl + X2 

^li[2*3*]^-4 

Xl - X2 



Xl +X3 



, - ~ +32 ?-l[3*2*]?-4 \ ,t t tj ,1, 
+ ?'l,3''2+4,4 44 d4a2a3di + h.c. 

X1-X3 J 



(A41) 



The symbols mean: tfj = XiJ"t''Xjc, and [a/3] = £q£^ + ^(£a ^ £/3)^^^- The star, *, means that the corresponding 



polarization vector is complex conjugated: * ^ s*. The quartic gluon term with derivative reads 

S[aAA]2 i234alaJ4a4 + h.c^ + X^qaa]^ I234ai4a3a4 

1 OQ/I 



1234" 



(A42) 
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where 



1 



-[aAA]2l234 = — g 

+ n+3,1^4,2 Ei^s • £2^4 

+ '~3+2, 3*^4, 1 £3^2 ' ^1^4 



n+2,ir4,3eie2 ■ £3^4 — — r~m / / 



{XV+X2Y 
{Xi - X3)(X2 + Xj) 

(X3 -X2)(X1 +3:4) ^^3^^ ^^^^^ 



(A43) 



[dAAY 1234 — ^ 



''l+2,l''3+4,3ei£2 ■ ^3^4 



{Xl + Xl)'^ 



- [?3,1?2,4 + ?1.3?4,2] e\e^ ■ ele, ^""^ j^^-^l)^ ""'^ /-i^s jae.c. 

- [r3,2ri,4 + r2,3r4,l] e\e^ ' ^2^4 + ^3)(xi + X4) _^„e.c._^ac.C3 

— a;4 j 

The instantaneous gluon interaction between quarks and gluons reads 



%AA](V.V) = E / [1234]^(pt - p) ^2 ^_^al2 ^a^ {} , 

1234 



(A44) 



(A45) 



where the brackets {} contain, 



{} = 

- ei£2 



+ eie2 



a;2 - Xl 

2 2(X1+X2)2 
Xl + X2 

(X3 + X4)2 
Xl + X2 



^^1+2,1^3+4,3 ^130-^X4 a\a\dzhi + fi+2,ir4,3 X3X4 63aia2^4 - ri+2,ir3,4 ^l3^-4 S^a\a\d^ 
^3+4,3^2,1 X3^^^-4 h\d\a\a2 
(^"2,1^4,3 + ^,2^3,4) ^13^-4 4«i«2d3 - (^2,1^3,4 + n,2r4,3) X3X4 h\a\aihi + /i.e. . (A46) 



(Xl - X2)^ 

Finally, the instantaneous gluon interaction between quarks reads 



^{■<1>^Y = Y^l I1234]5(pt 4^x1x2x3x4 {} , 

1234 



(A47) 



where the brackets {} contain, 
{} = - 



X1X2X3X4 ,a .a \~ ~ , ~ ~ 1 X3X2X1X4 ,a .a \~ ~ , - ~ 1 
2 ^12*34 ri,2r4,3 + ?'2,ir3,4j - -J^ ^32^14 r3,2r4,l + r2,3ri,4j 



(xi - x^f- 



(X3 - X2)2 



+ 



+ 



?-2?-l?-4?-3 .a .0. \~ ~ , ~ ~ 1 ?-2C-3?l4C-l .a .a \~ ~ , - ~ 
— *21*43 Fl,2r-4,3 + r-2,ir3,4j — ^ ^2 t23*41 [^,2^4,1 + r2,3ri,4. 



(xi - x^f 



(X3 - X2)2 



616362^4 

rfjrf3d2d4 



XlX2xk^?-4 .a .a X3X2xI '^'^^-4 .a-aj. ^ 

h2h4 '^2,ir'3+4,3 7 7^ *32*14 »^2,3?^1+4,1 



(Xl - X2)2 



(X3 - X2y 



^ *12*43 ^4,3r-l+2,l _|_ *13*42 '^4,2ri+3,l 



(Xi + X2)2 



b\bldlb2 + h.cj 
b\d2dld4 + h.c. 



- 2 



XlX2C-4^-3 ,a ,a r- - , - . 1 j,t j j. , o ^l^^^-3?-4^^X2 ,„ ,„ 
-*12*43ri,2»'4,3 + ''2,ir-3,4] 614^462 + 2 ■ 



(Xl - X2)^ 



(Xl + X3)^ 



*13*42 n+3,ir2+4,2 6l4c^4&2 • (A48) 



Useful color identities are: t^H'' = -t"/(2iVc), tH'' + fHH" = i{Nj2)t'' . 
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APPENDIX B: QaQa INTERACTION 



where 



Several factors are needed to estimate the small- a; be- 
havior of the potential kernel vq{13, 24) in Eq. H71|) . Mo- 
menta are labeled according to Figs, ^andl^l 



where 



with 



/l3,24 = exp [-(A^?3--^24)VA1 , 



M 



XiXi 



k^A — - hi^ 



(Bl) 



(B2) 



(B3) 
(B4) 



and. 



where 



\z\ 



Xi 



Xi 



Ml.-m^ = Vs., 



T. ^3 , . 



X3 



'■13 



(B9) 



(BIO) 



= (Bll) 

X3 



while the form factors /a, 54/51. 2 have the arguments, 



where 



X13 - -^24 = 4(fci3 + ^24) q, 



q = ki3- fc24 , 



(B5) 



(B6) 



is the momentum transfer that goes over to the standard 
one in the NR limit. 



(xi - Z){X3 + z) 



(B7) 



In the last term in Eq. (|71|) . the form factors /i, 52/53. 4 
have arguments 



Xl — z 



(B8) 



M^ 



X3 



54 



X3+ z 



■ V3 



(B12) 



and, 



Pi 



(B13) 



The last term in Eq. 17111 can be written with the coef- 
ficient 1 = (1 — /i3,24) + /i3,24 but only the second term 
counts at small z because (1 — /13.24) is proportional to 
the momentum transfer squared. The factor /i3,24 be- 
comes common to all terms in Eq. (I71f) and is taken out 
in front. The LF instantaneous term can be split into 
the part // that joins the low-energy exchange 1 — // 
that goes with the high-energy exchange. This way one 
obtains 



™2 ~2 



X1X3 



Xl X3 



M^U(k^3,Xi) 



^da;2d2^^ /^/l3^~^(^3^24)^(44,X2)=0, (B14) 



X2X4 xt, 



where 



■Do (13, 24) = 6'(^)^5/l''5/4 {v+Mgh + V+low) + d{- z)f^ /■g.fz /2 (U- Mgh + V-^low) ■ 



(B15) 



The gluon mass ansatz contributes to the low-energy exchange terms only. In terms of the invariant masses from Eqs. 
(lB8l-lBT3l. 



/i,52/53,4 -l \xl+xl {Ml2-m^){Mls-m^) 



/3,54/5i,2 Ul + xl {Ml,- m^)(Mli - m2) 



V+ high 



2 i X3X2 {Ml,-m^Y + {Ml^-w?Y 

These have the same limit when z — > for fixed q*^, 
fim v+ high = lim u_ high 



(B16) 
(B17) 
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/a, 54/51, 2 - 1 X1-X3 {q-^ 



^13 



(B18) 



The terms on the order of and higher are finite when divided by the square of X5 = \z\. Terms hnear in z produce 
an integral convergent in the sense of principal value 2, 16] . When ^ ^/z 0, (P'k,^^ removes one power of z from 
the denominator in Eq. (|7()|l . while v±high vanish for 2-^-0. The contributions of v±high are k"^ /m? times smaller 
than the dominant terms and can be ignored in the first approximation. One can see this by integrating v± high with 
a Coulomb wave function. 



The low-energy terms read 

/i. 52/53, 4 
V+low = -. 



m2)/a;4-/x2(2,5,3)/x5 



m2)/a;i -/i2(2,5,3)/x5 



52 



m^)|xl+^i■^{2,?>,i)/x^ 



53 



m^)|xi + ^i■^{2,^,i)/x^ 



(B19) 



/3, 54/51, 2 



{Ml, 



m2)/a;2 -Ai^(l,5,4)/x5 



54 



m2)/a;3 -Ai^(l,5,4)/a;5" 



{Ml, 



mP-)/x3 + ljfi{l,b,4)lx^ 



{Ml, 



TvP-)/X2+ ljfi{l,^A)/xs,_ 



(B20) 



APPENDIX C: MASS TERMS 

The mass terms with i = 1,3 in the eigenvalue Eq. 
(|ljl4|) are given in Eq. with given in Eq. (|^ . 

and A^§ in Eq. TOq originates from Eq. with 

A = Aq . Namely, the quark mass counterterm in X of Eq. 
(0) adds Srn^j to the original mass parameter in Eq. 
and the free ultraviolet-finite part of the counterterm 
is such that ml appears in Eq. 170(1 . The condition on 
ml that the eigenstates of H\g with quantum numbers 
of a single quark have eigenvalues growing to infinity is 
fulfilled below by representing gluon interactions in the 
case of the single quark state by a new gluon mass ansatz. 
The resulting value of TOq enters into the quarkonium dy- 
namics. The determination of the ultraviolet-finite part 
of the mass counterterm in X in Eq. is thus based on 
the picture that comes out from simultaneous consider- 
ation of two eigenvalue equations, one for the state with 
quantum numbers of a single quark (or an anti-quark, the 
result is the same), and another one for the quarkonium. 
The key physical assumption made in the comparison is 
that the binding of effective quarks in the quarkonium 
state occurs at the expense of change in their individual 
structure. While the buildup of self-interacting clouds 
of gluons around single quarks leads to the infinite quark 
masses, in the case of a colorless pair the main parts of the 
gluon clouds can recombine into a colorless object that 
may fiy out of the region of strong interaction with the 
quarks, leaving behind only the minimal remnants of the 
gluon clouds required to form the quarkonium eigenstate 
with a finite mass. The new finite balance is described 
using the gluon mass ansatz parameter 6^ . The finite bal- 
ance can be achieved because the quark-anti-quark state 
looks neutral from large distances and does not continue 
to generate gluons over infinite distances along the LF. 
This scenario is partly similar to the one originally de- 
veloped in the LF dynamics in [i^ l43l | , and studied in 



|43.l45| . The main differences are related to the fact that 
the physical picture that emerges here in the finite ef- 
fective theory with the gluon mass ansatz relies on the 
phenomenological parameter S^. A formal cutoff param- 
eter S'^ of the canonical theory, the coupling coherence 
phenomenon that may work over many scales of an ultra- 
violet cutoff, and the condition of transverse locality are 
not employed in the new picture. Instead, the present 
scenario can be studied in higher orders of perturbation 
theory according to the known rules [30I l33 | that explic- 
itly preserve the boost invariance, cluster decomposition 
property, and unitary connection with the initial theory. 

The eigenstate of Hxg with a single quark quantum 
numbers and momentum p with components p+ and , 
has the eigenvalue 

p- = {p^' + m^)/p+, (CI) 

and the decomposition in the effective particle basis, 

b) = IQA„) + |QA„<?A„)+ .... (C2) 

The new gluon mass ansatz is introduced in the quark- 
gluon component. It is different than in the quarkonium 
case because the states have different quantum numbers 
and are made of different numbers of effective particles. 
Dropping the subscript Aq as in Eq. ((Sl|l . the eigenvalue 
problem is written as 

iT, + fg)\Qg) + Y\Q) = E\Qg), 

Y\Qg) -f T,\Q) = E\Q). (C3) 

The new ansatz enters through the kinetic energy Tg, 
which contains 

where x and refer to the relative motion of the ef- 
fective gluon with respect to the quark. The operator 
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R from Eq. (|57|l is now replaced by the one with P 
that projects on the single effective quark basis state with 
kinematical momentum components p"*" and p^. In the 
perturbative expansion in g, only the second term on the 
right-hand side of Eq. ljC4p contributes in orders up to 
g^. Thanks to the boost invariance, the resulting eigen- 
value condition reduces to an equation for m^, which is 
independent of p, 



37r2 



M'q — rri^ 



For fh? to be positive and growing to infinity when (5^0, 
one can write mg in the integral form, 



1^ jdxdPn^fl{x) 



1 A^2_^2 



X"^ Alg — TV? 



with some function A^H that satisfies the condition, 



M'^n > Ml >m^ 



This condition can be satisfied by writing. 



Mi = K ' + ^ii{x, K^)]/x +{k^' + m')/{l - x) 



(C8) 



(C9) 



(C5) and assuming that 



Mq > Mo > . 



(CIO) 



^ ^ As long as the difference /ig — /^.g does not vanish for 
X — > 0, the single quark mass will tend to oo when (5 — s- 0. 
But this may easily happen here because the larger is the 
gluon mass ansatz /ig , the stronger the single quark mass 
eigenvalue diverges in the limit 5 ^ Q, while /ig remains 
free to vanish in the limit a; — > and lead to a finite 
mass contribution in the quarkonium dynamics. Using 
(C7) Eq. (jnH) for ml one obtains Eq. ((EHJ- 
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